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Abstract. An important question about resonance extraction is how much resonance poles 
and residues extracted from data depend on a model used for the extraction, and on the precision 
of data. We address this question with the dynamical coupled-channel (DCC) model developed 
in Excited Baryon Analysis Center (EBAC) at JLab. We focus on the Pn ttN scattering. We 
examine the model-dependence of the poles by varying parameters to a large extent within the 
EBAC-DCC model. We find that two poles associated with the Roper resonance are fairly 
stable against the variation. We also develop a model with a bare nucleon, thereby examining 
the stability of the Roper poles against different analytic structure of the Pn amplitude below 
ttN threshold. We again find a good stability of the Roper poles. 



1. Introduction 

Extraction of N* information, such as pole positions and vertex form factors, is an important task 
in hadron physics. This is because they are necessary information to address a question whether 
we can understand baryon resonances within QCD. In order to extract the N* information, first, 
one needs to construct a reaction model through a comprehensive analysis of data. Then, pole 
positions and vertex form factors are extracted from the model with the use of the analytic 
continuation. Therefore, the N* information extracted in this manner is inevitably model- 
dependent. There are several different approaches to extract the N* information. Although 
almost all 4-stars nucleon resonances listed by Particle Data Group (PDG) are found in all 
approaches, existence of some N* states, in particular those in the higher mass region, is 
controversial. Thus, commonly asked questions are how much model-dependent the extracted 
resonance parameters are, and how precise data have to be for a stable resonance extraction. 
These are the questions we would like to address in this workpp, within a dynamical coupled- 
channels model (EBAC-DCC) \2\. We focus on the irN Pn partial wave and the stability of 
its pole positions, particularly those corresponding to the Roper resonance. In the region near 
Roper N (1440), two poles close to the irA threshold were found in our recent extraction [3J from 
the model obtained by the fit to irN — > irN scattering data (JLMS) [3]. 

while only one pole in the similar energy region was reported in some other analyses. We 
examine the stability of this two-pole structure against the following variation, keeping a good 
reproduction of SAID single-energy (SAID-SES) solution [5J unless otherwise stated: 



• Large variation of the parameters of the meson-exchange mechanisms as well as bare N* 
parameters of the EBAC-DCC model. 

• Inclusion of a bare nucleon state: The analytic structure of this model is rather different 
from the original EBAC-DCC model, in particular in the region near the nucleon pole [6], . 

• Fit to the solution based on the Carnegie-Mellon University-Berkeley model (CMB) [7] 
which has rather different behavior from SAID-SES for higher W. 



2. Dynamical coupled-channels models and analytic continuation 

In this section, we briefly describe dynamical coupled-channels model used in this work, followed 
by a brief explanation for the analytic continuation used to extract poles from the model. 



2.1. EBAC-DCC model 

The EBAC-DCC model contains irN, rjN and irirN channels and the irirN channel has irA, 
pN and o~N components. These meson-baryon (MB) channels are connected with each other 
by meson-baryon interactions (vmb,M'B'), or excited to bare N* states by vertex interactions 
(Tmb-^n* ) • With these interactions, the partial-wave amplitude for the M(k) + B(—k) — > 
M'(k') + B'(-k') reaction can be written by the following form: 

TMB,M'B'(k,k ,E) = tMB,M>B> (k,k',E)+t MB 

,M'B' 

(k,k',E), (1) 

where the first term is obtained by solving the following coupled-channels Lippmann-Schwinger 
equation: 

tMB,M'B'(k,k',E) = V M B,M'B>(k,k') 

+ / q 2 dqv M B,M"B"(k,q)GM"B"(q,E)t M "B",M>B>(q,k',E).(2) 

M"B" M" B" 

Here Cmb is the integration contour in the complex-q 1 plane used for the channel MB. The 
second term of Eq. ([1]) is associated with the bare N* states, and given by 

^MB,M'B'{k: 

k',E) = £ 

Fmb-±n* (k, E)[D(E)]ijT N *^ M > B >(k' , E), (3) 

hi 

where the dressed vertex function T n*^m' B'{k, E) is calculated by convoluting the bare vertex 
^N*^,M'B'{k) with the amplitudes tMB,M'B'(k,k' , E). The inverse of the propagator of dressed 
N* states in Eq. © is 

[D-\E)} i!:j = (E — m N *)5 itj - Xi,j(E), (4) 
where Tin is the bare mass of the i-th TV* state, and the TV* self-energy is defined by 

i 

Sij(E) = ^2 / q 2 dqf N *^MB(q,E)GMB(q,E)r M B->N*(q,E). (5) 



MB 



Defining E a (k) = [m^, + k 2 } 1 / 2 with m a being the mass of particle a, the meson-baryon 
propagators in the above equations are: GMB^k, E) = 1/[E — E]\i(k) — Es(k) + ie] for the 
stable ttN and r]N channels, and Ga/b(A;, E) = 1/[E — EM(k) — Es(k) — Y<MB(k, E)] for the 
unstable ttA, pN, and o~N channels. The self energy Y>MB(k, E) is calculated from a vertex 
function defining the decay of the considered unstable particle in the presence of a spectator tt 
or N with momentum k. For example, we have for the irA state, 

y (h _ ™, A f 2 , M nN (q) I/a^aK 1 ?)! 2 / r n 

**Al«, *j EA(k) q [M 2 N (q) + ^]V2 E _ E ^ k) _ [A ^ (g) + fc2] i/2 +ie > W 

where M 7r jv(Q) = E n (q) + E]y(q) and ./a-s.ttaK'?) defines the decay of the A — > irN in the rest 
frame of A, C3 is the corresponding integration contour in the complex-g plane. 



2.2. Bare nucleon model 

To examine further the model dependence of resonance extractions, it is useful to also perform 
analysis using models with a bare nucleon, as developed in Ref. |8j. Within the formulation 
given in Sec. 12.11 such a model can be obtained by adding a bare nucleon (Nq) state with mass 
mjy and Nq — > MB vertices and removing the direct MB — > N — >■ M'B' in the meson-baryon 
interactions vmb,M'B' • All numerical procedures for this model are identical to that used for the 
EBAC-DCC model, except that the resulting amplitude must satisfy the nucleon pole condition: 

t*N,irN(k -»■ ^on, k -> k on , E -> m N ) = - ~- -. (7) 

hi — m N — 2j[m]\[) 



with 



m N = m ° N + t(m N ) and F 7rNN (k on ) = F^(fc on ) . (8) 



Here we have used the on-shell momentum defined by E = J m? N + k% n + \/ m\ + A^ n (Im[A; on ] > 

0). Also, S(mAr) is the self-energy for the nucleon. More details for the calculational procedure 
following Afnan and Pearce are found in Refs. [HE]- 



2.3. Analytic continuation 

Once a fit is obtained, we then apply the method of analytic continuation to find resonance 
poles. The procedures for performing this numerical task have been discussed in Ref. [31 [9]. 
To search for resonance poles, as discussed in the above references, the contours Cmb and 
C3 must be chosen appropriately to solve Eqs. jS])-© for E on the various possible sheets of 
the Riemann surface. We only look for poles which are close to the physical region and have 
effects on the irN scattering observables. All of these poles are on the unphysical sheet of the 
ttN channel, but could be on either unphysical (it) or physical (p) sheets of other channels 
considered in this analysis. We will indicate the sheets where the identified poles are located by 
(sttTV, s v ]y, s^Ar, s-ttA, SpN, s a N), where smb and s nn iy can be u or p. The errors of the resonance 
parameters are estimated by using all values obtained in all fits we have performed. 



3. Results 

Now we show our numerical results to examine the stability of the P\\ poles. First of all, we show 
P\i amplitudes from JLMS and SAID-EDS (energy-dependent) [5] compared with SAID-SES in 
Fig. [TJ In Table [lj the pole positions from JLMS and SAID-EDS as well as x 2 P er data point 
(Xpd) are given. In the following subsections, we present results from various fits by varying 
the dynamical content of the EBAC-DCC model, using a model with a bare nucleon, and using 
different empirical amplitude for the fit. 



3.1. 2N*-3p and 2N*-4p fits 

We varied both the parameters for the meson-baryon interactions (vmb,m'B') and parameters 
associated with bare A^* states (m^», Tiyt^MB)- The obtained meson-baryon interactions are 
quite different from those of JLMS. We obtained several fits which are different in how the 
oscillatory behavior of SAID-SES amplitude for higher W is fitted. The results from the 2A^*-3p 
(dotted curves) and 2A*-4p (dashed curves) fits are compared with the JLMS fit (solid curves) 
in Fig. [2j The resulting resonance poles are listed in the 3th and 4th rows of Table [TJ Here we 
see again the first two poles near the 7rA threshold from both fits agree well with the JLMS fit. 
This seems to further support the conjecture that these two poles are mainly sensitive to the 
data below W ~ 1.5 GeV where the SAID-SES has rather small errors. However, the 2A*-4p fit 
has one more pole at Mr = 1630 — i45 MeV. This is perhaps related to its oscillating structure 



Table 1. The resonance pole positions Mr for P\\ [listed as (ReMjj,, — ImMp,) in the unit 
of MeV] extracted from various parameter sets. The location of the pole is specified by, e.g., 
(sjtN, s v n, s n7T N, s n A, s p n, s^n) = (upuupp), where p and u denote the physical and unphysical 
sheets for a given reaction channel, respectively, Xpd ^ s X 2 P er data point. 
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SAID-EDS 
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2iV*-3p 


(1368, 82) 
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2iV*-4p 
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2iV*-4p-CMB 
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Figure 1. The real (left) and imaginary (right) parts of the on-shell P\\ amplitudes as a function 
of the ttN invariant mass W (MeV). The solid curves are from the JLMS fit; the open circles 
are the SAID-EDS [5]. A is unitless in the convention of Ref. [5]. 



near W ~ 1.6 GeV (dashed curves), as shown in the Figs. [2j On the other hand, this resonance 
pole could be fictitious since the fit 2iV*-3p (dotted curve) with only three poles are equally 
acceptable within the fluctuating experimental errors. Our result suggests that it is important 
to have more accurate data in the high W region for a high precision resonance extraction. 

3.2. !N lN*-3p 

Here we show our results obtained with the bare nucleon model, and then address the question 
whether difference in the analytic structure of the TriV amplitude below TriV threshold strongly 
affects the resonance extractions. The bare nucleon model is fitted to SAID-SES, and at the 
same time, to the nucleon pole conditions Eq. (JSj) - Meanwhile, the original EBAC-DCC model 
has different singular structure below the ttN threshold. The question is whether such differences 
can lead to very different resonance poles. Our fit of the bare nucleon model is shown in Fig. [3] 
and compared with SAID-SES and JLMS. We see that the two fits agree very well below W = 1.5 
GeV, while their differences are significant in the high W region. The corresponding resonance 
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Figure 2. The real (left panel) and imaginary (right panel) parts of the P\\ amplitudes. 
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Figure 3. The real (left panel) and imaginary (right panel) parts of the P\\ amplitudes. 



poles are given in Table [TJ We also see here that the first two poles near the ttA threshold are 
close to those of JLMS. Our results seem to indicate that these two poles are rather insensitive 
to the analytic structure of the amplitude in the region below irN threshold, and are mainly 
determined by the data in the region m^r + m n <W< 1.6 GeV. 

3.3. 2N*-4p-CMB fit 

To further explore the dependence of the resonance poles on the data, we consider a solution 
from CMB collaboration [7]. This solution differs significantly from the SAID-SES mainly at 
W > 1.55 GeV. For our present purpose of investigating the stability of the lowest two poles 
near the ttA threshold, we fit the data which is obtained from replacing SAID-SES in the high 
W > 1.55 GeV region by the CMB solution. The results (dashed curves) from this fit are 
compared with JLMS in Fig. [TJ We see that the CMB solution has oscillating behavior near 
W ~ 1.6 GeV and this could be the reason why the fit has an addition pole near W ~ 1.6 GeV, 
as seen in 5th row of Table [TJ The large differences from JLMS at high W make the poles near 
W ~ 1.9 GeV very different; in particular their imaginary parts. On the other hand, their lowest 
two poles near the ttA threshold are close to other fits discussed so far. This again supports the 
above observation that these two poles are determined only by the data below W < 1.5 GeV 
which are reproduced very well in all fits. 




4. Conclusion 

We have examined the stability of the two-pole structure of the Roper resonance. We showed 
that two resonance poles near the n A threshold are stable against large variations of parameters 
of meson-exchange mechanisms within EBAC-DCC model [2]. This two-pole structure is also 
obtained in an analysis based on a model with the bare nucleon state. Our results indicate that 
the extraction of P\\ resonances is insensitive to the analytic structure of the amplitude in the 
region below irN threshold. We have also fitted to the old CMB amplitude, which is rather 
different from SAID-SES for W > 1.5 GeV, and still found that the Roper two poles are stable. 
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